Quantum criticality with multiple dynamics 
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Quantum critical systems with multiple dynamics possess not only one but several time scales, 
Ti ~ £ Zi , which diverge with the correlation length £. We investigate how scaling predictions 
are modified for the simplest case of multiple dynamics characterized by two dynamical critical 
exponents, z > and z<. We argue that one should distinguish the case of coupled and decoupled 
multiple dynamic scaling depending on whether there exists a scaling exponent which depends on 
both Zi or not. As an example, we study generalized 3> 4 -theories with multiple dynamics below 
their upper critical dimension, d + z < < 4. We identify under which condition coupled scaling is 
generated. In this case the interaction of quantum and classical fluctuations leads to an emergent 
dynamical exponent, z e = ^— j" z< ^ +1 ■ 

PACS numbers: 71.10.Hf, 71.27. +a, 73.43.Nq 



Quantum criticality has become a ubiquitous theme 
in condensed matter physics. The anomalous thermo- 
dynamics and transport properties close to a quantum 
critical point are fascinating as they defy paradigmatic 
theories of normal metals and insulators, and they are 
now frequently invoked to explain the unusual behavior 
observed in a wide class of systems ranging from heavy- 
fermion compounds, organic materials and high-T c su- 
perconductors to ultracold atoms. 

Quantum critical points are not only characterized by 
the symmetries of the system and the critical degrees of 
freedom, but also by their specific dynamics 1,2 . How- 
ever, often these attributes are not easily identified and 
simplifying assumptions are applied. For example, the 
quantum critical scaling hypothesis, that is commonly 
employed for the interpretation of experiments, assumes 
the critical behavior to be dominated by a single time 
scale r that diverges together with the correlation length, 
£, according to r ~ £ z , where z is the dynamical expo- 
nent. At finite temperatures T, scaling then implies the 
presence of a thermal length scale, £t ~ T~ x / Z , char- 
acterizing the distance beyond which thermal dominate 
over quantum fluctuations. However, often various crit- 
ical degrees of freedom with their respective dynamics 
coexist giving rise to multiple thermal lengths. In partic- 
ular, in quantum critical metals critical collective modes 
interact with critical fermionic quasi-particles that are 
both governed by different dynamics. In such cases, scal- 
ing assumptions become ambiguous as certain physical 
observables could be determined by one or the other dy- 
namics of the different critical degrees of freedom 2 - - — . 

Here, we consider quantum phase transitions where 
the order parameter is characterized by two dynamical 
exponents z > and z <1 with z> > z<, and study the 
scaling behavior of thermodynamics. For renormaliza- 
tion group (RG) fixed-points with multiple z, the free 
energy exhibits only weak or multiple dynamic scaling. A 
similar situation arises in the context of classical criti- 
cal phenomena^—, for example, for model C dynamics 
within the Hohenberg-Halperin classification^ where the 
order parameter couples to a conserved density field. The 



dynamical exponents together with the spatial dimension 
d give rise to two effective dimensions, D > = d + z> and 
D < = d+z < . Dimensional analysis then requires that the 
critical part of the free energy density, J-~ cr , is necessarily 
characterized by two scaling functions, /> and /<. The 
simplest generalization of the scaling hypothesis for F CI 
in the presence of multiple dynamic scaling is then given 
by (in the absence of dangerously irrelevant operators) 

7er(r,T) = J->(r,r)+J-<(r,r) 

= b- {d+z > ) f > {rb\ / \t < b\<,t > bl>) (1) 



b-^Mrb 1 ^,^, 
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with arbitrary scaling parameters b\ and bi- The scaling 
functions, /> and /<, depend on the tuning parameter 
r, that measures the distance to the quantum critical 
point and scales with the correlation length exponent v, 
and, generally, on two temperature scaling fields U = 
r)iT, with i G {>, <}, where T is the temperature and r\% 
are "kinetic coefficients". There are three characteristic 
length scales, £ ~ 61,2, where each of the three arguments 
of the scaling functions reaches a value of order one, 



j.-l/z< 



(2) 



The first can be identified with the correlation length at 
zero temperature, and the remaining are the two thermal 
lengths with £^ < £$. 

In the limit of small temperatures, T — > 0, but finite 
r, we can choose the scaling parameters such that b\ — 
62 = \r\ ~ v ~ £ r and the free energy simplifies to 

J- cr (r,T) = |rr(* f *>)/ > («7 ) t>|r|- , "> ) t<|r|-'*<) 

+ M" ( * f *<>/<(<7, t> |r|-"« > , t< |r|- w *< ) (3) 

with <j = sign r. At T = 0, the most singular part is then 
attributed to the contribution with the smaller dynamical 
exponent z < , 



7- cr (r,0)~|rr^<)/<((7,0,0). 



(4) 



So we arrive at the conclusion that the susceptibility 
Xrr = —d^J^ci) that determines, e.g., the correction to 
the compressibility for pressure tuned quantum critical 
points, is at T = determined by the mode with the 
smaller effective dimension d + z < . 

At finite temperatures T > 0, the generalized scaling 
hypothesis (JT|), however, is only of limited use in predict- 
ing the critical behavior, in particular, in the regimes 
6- > £t- Consider the free energy at scales of the 



two respective thermal lengths, b\ = t 
bo=t. 



-l/z< 



= (*>)^7> (rt^,t<(t>)-^,l 



£p and 



(5) 



d + z< 



7< 



) «<,l,i > (t<)- 



The temperature dependence is here not explicit due to 
the residual dependence of the two functions ft on tem- 
perature, so that the critical behavior is not evident. In 
particular, the third argument of the /< function diverges 
for T — > 0, and the critical behavior (r = 0) of thermo- 
dynamics will crucially depend on the analytic properties 
of/<- 

We can distinguish two cases. In the case of decoupled 
multiple dynamic scaling the asymptotic behavior of the 
two functions /< and /> in Eq. ([T]) does not dependent 
on the respective other thermal fields, i> and t <: so that 
the temperature dependence becomes manifest and scal- 
ing laws follow with exponents involving either z > or z < 
depending on the quantity of interest. In the case of cou- 
pled multiple dynamic scaling, however, the dynamics of 
the two modes mix so that the asymptotic behavior of the 
two scaling functions does depend on both, t> and i<, 
and Eq. is not sufficient to predict the temperature 
dependence of thermodynamics. In particular, scaling 
laws can then arise with unusual exponents depending 
on both dynamical exponents z > and z< . 

In order to identify a mechanism leading to coupled 
multiple dynamic scaling, we examine a specific exam- 
ple in the following. We consider the effective theory 
for a d-wave Pomeranchuk (or nematic) instability in 
two-dimensional isotropic metals describing a sponta- 
neous deformation of the Fermi surface. The order pa- 
rameter is a quadrupolar tensor with two polarizations, 
<1> T = (</>i, 02), effectively described by a $ 4 -theory 



S 



d d xdr 
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The dynamics is generated by the coupling of the col- 
lective quadrupolar mode to particle-hole pairs in the 
metal, and, interestingly, it differs for the two polariza- 
tion of the order parameter. The quadrupolar fluctua- 
tions at a given momentum k are damped by particle- 
hole pairs in the metal if the polarization is longitudinal 
to kikj — Sijk 2 /2 leading to a dynamic exponent z > = 3 
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FIG. 1: One-loop corrections to (a) the tuning parameter and 
(b) the interaction u. 



while it remains undamped for the transverse polariza- 
tion giving z < = 2--. While the renormalization group 
flow of this model is towards a stable Gaussian fixed- 
point, it was demonstrated in Ref. [6] that the multi- 
ple dynamics results in an extended quantum-to-classical 
crossover. Moreover, it was observed that this extended 
crossover leads to logarithmic corrections to thermody- 
namics reminiscent of the coupled multiple scaling dis- 
cussed above. In order to study the coupling of dynamics 
in the presence of an interacting fixed-point we generalize 
the model to an effective dimension D — 4 — e which sta- 
bilizes the RG flow towards a Wilson-Fisher (WF) fixed 
pointii so that universal physics at lowest energy emerges 
and the properties of the scaling functions can be studied. 

The Green function in ^ possesses the following mo- 
mentum and frequency dependence^ 



G 1 (&u n ) 









ff> X (<f,W n ) 



(7) 



where the eigenvectors itself depend on the orientation 
of momentum q via the rotation matrix 



cos (20) 



.(20)\ 



sin (26) cos (26) J ' 



(8) 



that depends on the angle between momentum and, say, 
the x-axis, 9/L(q, x), and reflects the d-wave symmetry of 
the order parameter. In the disordered phase, ($} = 0, 
the respective Green functions are given by 



9<(q^n) = r + q 2 +7f < - 



,2z°-2 ' 



(9) 



-It- \ , 2 , 

9> (q,u n )=r + q + ??> y0 _ 2 , 



with q = \q\. The tuning parameter of the quantum 
phase transition is r, and rji, with i £ {>, <}, are kinetic 



coefficients. We consider effective dimensions d ■ 
and d + z< 



> 4 



4 — e so that the whole theory is below its 
upper critical dimension, = 4, due to the z° -mode. 



We will use z2. = 3 and discuss the two realizations d = 2, 
e and d 



z 



ii 



< 



e, z% 



2. The original model is 
recovered in the limit e — > 0. 

We perform a one-loop RG analysis by subsequently 
integrating out high energy modes within a momentum 
shell (A/6, A) with logfr > 1 and A being a momentum- 
cutoff that we set to unity in the following, A = f. 
The form of the theory is restored by rescaling momenta 
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q — > q/b and frequencies w„ — > Lo n /b z . To ensure an unbi- 
ased treatment of the multiple dynamics, we use here an 
unspecified dynamical exponent z. We will later demon- 
strate that physical observables do not depend on the 
specific choice of z. As a result of this scaling scheme, 
the kinetic coefficients, rji, obey in one- loop order the 
tree-level scaling equations = (zf — z)T]i, where 

i E {>,<}. The quantum fluctuations, see Fig. [TJ lead 
to an RG flow of the coupling u and the tuning parameter 
r, which at temperature T — is given by 



Or 



d\ogb 

du 
d log b 



2r 



1 



24ttt] < 
(4-d-z)u- 



327r?7< 



(10a) 
(10b) 



The appearance of the kinetic coefficient rj < indicates 
that the RG flow at T = is determined by the loop 
corrections attributed to the z° -mode that possesses the 
smaller effective dimension d + z< = 4 — e, while the 
z > mode gives only subleading corrections. Introducing 
the coupling U = u/rj < , one finds that for scales b > 
&wf = EAvf/^o — l| 1//e w ith the bare coupling Uq the 
interaction is driven towards the WF fixed-point U —> 
Uwf = 327re/3. Furthermore, from the scaling equation 
for the tuning parameter follows the correlation length 
exponent 



1 

v 



(11) 



Note that this value for the correlation length exponent 
does not belong to the standard universality class of 
O(N) models with diagonal propagators 14 . In the fol- 
lowing, we consider the theory at scales b 6wf where 
universality emerges and study the RG flow in the vicin- 
ity of the WF fixed-point. 

Temperature T is a relevant scaling field with the scal- 
ing dimension given by z. For our purposes it is actually 
convenient to consider the scaling of the two temperature 
fields t y = r/yT and t< = ?7<T, 



dU 



d log b 



Zy U 



and 



d log b 



z<- 1< 



(12) 



with the identification 



at one-loop order. The 



scale where either of the two fields reaches the cutoff 
scale, ti(bj.) = 1, i 6 {<,>}, defines the two ther- 
mal lengths, b l T — introduced in Eq. @. At the 
thermal length scale, the so-called quantum-to-classical 
crossover— occurs where critical modes change their char- 
acter. Whereas for small scales, i.e., for large momenta 
q > ^y 1 , the critical modes have an essential quan- 
tum character with effective dimension d + z, for scales 
q^ 1 > £t they become classical with effective dimen- 
sion d only. Due to the presence of multiple thermal 
lengths, however, the quantum-to-classical crossover in 
the present case is extended^. For scales b between the 
two thermal lengths, b^ < b < 6^, the quantum fluctu- 
ations of the z<-mode coexist with the thermal fluctua- 
tions of the z>-mode, see Fig. [51 Interestingly, both quan- 
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FIG. 2: Quantum-to-classical crossovers of the 2>-mode and 
z<-mode as a function of momentum. There exists an ex- 
tended quantum-to-classical crossover regime for £p> > q > 
£t< w here quantum fluctuations of the z<-mode coexist with 
thermal fluctuations of the z>-mode (gray shaded area). 



turn and classical fluctuations contribute to the RG flow 
within this range, see Fig. [TJ so that for b > b^, ^> 6wf 



dr 
dlogb 



2r 



247Tf7< 



uT 
6ir 



1 

— r - 
v 



tAvF t< 

(in 



(13) 



The quantum fluctuations of the z<-mode determine v 
as before. The classical fluctuations of the z>-mode lead 
to an additional term linear in T. This term is cru- 
cial as it induces a new length scale that reflects the 
interplay of the different dynamics. The coexistence 
of quantum and classical fluctuations mixes the scal- 
ing fields r and t < giving rise to a new scaling field 

R 0>) = K&) + 6 / { iZ%) *<(&) which obe y s 



d log b v 



(14) 



The condition R(£r) = 1 allows to identify the new 
length scale, £r = (i?(6^)) _!y 6^, that reads explicitly 



16e 



v r\ < (1- 



aT 1 - 



i/z<)/(i/z>) 



with a — ,. . 

9 1 — z< 

'emerging' dynamical exponent given by 



(15) 



This defines a new 



u(z> - z<) + 1 



(16) 



This simple formula encodes in a nutshell how the inter- 
play of the two time scales can lead to qualitatively new 
physics. At first glance it is surprising that all explicit 
dependence on the number of space dimensions cancels 
despite the fact that the contribution arises from classical 
fluctuations. The information on the number of dimen- 
sions is encoded in the RG equation for U. While uT 
has the engineering dimension 4 — d, close to the Wilson- 
Fisher fixed-point it gets replaced by UwFt< with dimen- 
sion z < . For the same reason, there is also no explicit 
dependence on the interaction strength despite the fact 
that d + Zy > 4. 

The generated temperature dependence in Eq. (|15p re- 
sults, in principle, in the presence of three length scales, 
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d = 2, 4 = 2 - e 


d = 2 - e, z< = 2 


z e 








4f e 




Xtt ~ 


rpd/ Z> — 1 rp — 1/3 


j>d/z>-l _ j.-(l+e)/3 


XrT ~ 


rp(d— 1/v+z<:)/z B -l _ rp7e/Si 


y(d-l/y)/z< _ jr-5e/18 


Xrr ~ 


rp(d+z < —2/i/)/z e _ rp — e/18 


j.(d+z<-2/i/)/z< _ m-e/18 



TABLE I: Emergent dynamical exponent z e , see Eq. (|16p . for 
two choices of d and z<, and the generalized susceptibilities 
at criticality r — as a function of temperature T. Coupled 
multiple dynamic scaling obtains for z e > z< with scaling 
laws depending on z e . 
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FIG. 3: Quantum critical phase diagram and crossover lines 
(dotted) in the (r, T) plane for (a) decoupled and (6) coupled 
multiple dynamic scaling. The shaded area is bounded by 
the Ginzburg line. While the specific heat is sensitive only to 
the lower and the compressibility to the upper crossover line, 
the thermal expansion changes its behavior at both crossovers 
(see appendix). 



£> ~ T~V Z >,£< ~ T-!/ z < and £ R , even at critical- 
ity r = 0. The generated length £r ~ T -1 / 2 " obeys 
£r/£t > 1 as always z e < z> but the ratio £a/£t de- 
pends on the relation between the emerging exponent z e 
and z < . Its asymptotic behavior distinguishes whether 
the multiple dynamic scaling turns out to be coupled or 
decoupled. This is best seen by considering the free en- 
ergy density of Eq. (QJ . For z e < z < the RG flow at r = 
is cut off by temperature, T, before £r is reached. As a 
result the free energies Ti ~ {^ r )~ l - d+Zi \ i 6 {<, >} are 
determined by their respective thermal length, £^ and 
so that the multiple dynamic scaling is decoupled, 
see appendix for details. In contrast, for 

z e > z < 1 > vz < <^ 1 > vz e (17) 

realized for dimensions d > 2 — |e, the two length 
scales £r and £^ govern the critical properties so that 
J-> - (^>)-( d + z >) and J"< - (£ R y (d+z<) - In particu- 
lar, the part T< is determined by and thus depends 
on the emergent dynamical exponent z e resulting in cou- 
pled multiple dynamic scaling. Table U collects the results 
for the generalized susceptibilities \ab = —d a dbJ- OT with 
a, b = T, r at criticality r — 0. For pressure-tuned quan- 
tum criticality, XrT and \rr can be identified with the 
thermal expansion a and the critical part of the com- 
pressibility, k, respectively, and xtt is just the specific 
heat coefficient, 7. If z e > z<, we obtain coupled multi- 
ple dynamic scaling (i.e. scaling exponents which depend 
on both z < and z>) for the thermal expansion, XrT, and 
the compressibility, Xrr- 

Comparing the three thermal length scales with the 
correlation length at T = 0, £ r ~ identifies vari- 

ous crossover scales in the phase diagram for r > 0, see 
Fig. [31 where the thermodynamics changes its behavior 
(for more details, see appendix). For coupled multiple 
dynamic scaling, one crossover line depends on z e . As 
usual, for a sufficiently negative r an additional crossover 
occurs from the quantum WF fixed-point to a classical 
one giving rise to the Ginzburg line in Fig. ([3]). As the 
model ((6]) reduces in the classical limit to a XY-modcl the 
classical phase transition is of Kosterlitz-Thouless type 
for d = 2. 



To summarize, we studied quantum critical thermody- 
namics in the presence of multiple dynamical scales, and 
we distinguished between decoupled and coupled multi- 
ple dynamic scaling scenarios. Coupled scaling occurs 
when the classical fluctuations arising from the mode 
with the larger dynamical critical exponent generates the 
mass term which then grows rapidly due to quantum fluc- 
tuations of the other mode. This quantum 'boost' leads 
for vz < < 1 to an emerging thermal length T~ x / Z ' with 
the exponent z e of Eq. (flfl)) . Note that this exponent and 
the mechanism how it is generated are - at least in one- 
loop order - completely universal. z e depends only on 
the two dynamical critical exponents z < and z> and on 
the correlation length exponent v but not, for example, 
on the number of space dimensions. As it is also inde- 
pendent of the specific choice of propagators of Eq. ^ , 
we believe that it is relevant for a wide class of models. 
The additional length scale T _1 / Ze gives rise to unusual 
scaling-laws in thermodynamics. To search for it exper- 
imentally, one option is to look for phase diagrams as 
shown in Fig. 13^6) where an upper crossover line is con- 
vex (vz e < 1) while a second crossover line has a larger 
exponent. 
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Appendix 

The free energy density J- per volume V is defined as 
usual via the temperature T and the partition function 
Z as T = —(T/V) InZ. The partition function is deter- 
mined by the functional integral Z = J P$ exp(— S) in 
terms of the action S defined in Eq. (6) of the main text. 

Introducing a momentum cutoff \q\ < A, we can derive 
an RG trajectory for the free energy similarly to Refs. [HI, 
1 1 61 ] . Integrating out a high energy shell \q\ G [A/6, A] with 



5 



ln(6) <C 1 yields a contribution to the free energy density J r > , with ti = r\i T, we obtain 



^ J|9|e[A/6,A] 



In P(r + q z + 



2 , (V<^n) 2 



n d A d b-i 

(2TT) d b 



duj , f u) \ 
mth ( 1 



+ In [ /3(r + q + z0 ^ 



Tin 2sinh 



(A. la) 



2tt 



coth ( — J arctan 




(A.lb) 



with fid being the <i-dimensional solid angle. A differen- 
tial RG equation for T follows by considering the limit 
b — > 1 + . Furthermore, we note the free energy density 
has the engineering dimension d + z, with z being the 
arbitrary dynamical exponent used for the scaling of the 
temperature, di n ^ T — zT, and thus T(b) — Tb z . Solv- 
ing the RG equation one obtains the total free energy 
density in the form of an integral over an RG trajectory. 
In one loop order, the latter reads 



dU ^^{b) 



(A.2) 



where 



dbb- 1 b-^ d+z ^ 



A d n d r°° du 



x arctan 



(2n) d J Q 2n 



( u \ 



A z >- 2 (r(6) + A 2 ) 



coth I — — I 
\2T(b)J 

(A.3a) 



fOO A d (") 

F <=1 dbb-H-^^T(b) 



x In 2 sinh 



(27r)° 

y/A 2z <- 2 (r(b) + A 2 )" 
2rj<(b)T(b) 



(A.3b) 



As discussed in the main text, the free energy den- 
sity T has two contributions stemming from the two 
modes. The scale dependence of the temperature T(b), 
the kinetic coefficients f]<,>(b) and the mass r(b) were 
defined in the main text. In the following, we denote 
e.g. by T(b) the running temperature and by T without 
the explicit b dependence the starting value for the RG 
flow. Redefining first u — eT(b) and then substituting 

A/6 = t x l z< e~ l for the part T< and A/6 = t% z> er l for 



,(d+z > )/z > r co (-> 

J"> = ^ / dle~ dl [td 




(2ir) d ./n 2?r 



de /e 
— coth - 
2 



c e 



r(A<> 1/2> e')/A 2 + 1 



(d+ 2< )/4 




(A.4a) 



(A.4b) 



dle- dl -^-j 

\Jr{At~ 1/z< e l )/A 2 + l 
2e z < 1 



Note that with these substitutions all dependencies 
on the arbitrarily chosen scaling exponent z have van- 
ished. Next, we write Eqs. (IA.4I) in a scaling form. 
For T > r vz> , the flowing mass can be expressed in 
terms of the emergent scaling field R(b) and the constant 
a = {Wwf)/(67t(1 _ z < i/)), 



r(b) = R(b) -a v< (b)T(b) 



(A.5) 



With R(b) = Rb x l v = A 2 t B 1/v {A/b)- 1 /" and 
r] < (b)T(b) = t < b z < (see Eqs. $T2$) and ((T4j) of the main 
text) we find that the free energy density J- < depends on 

A- 2 r(Ai< 1/z< e l ) (A.6) 

and therefore only on the combination t < l ^ uz 
can thus be expressed in the scaling form 



R. It 



T K = 6-( d+z <) — f < (Rb 1 ^,t < b z <) 



(A.7) 



From the usual scaling analysis one can now deduce that 
the free energy density of the mode with the smaller dy- 
namical exponent scales for T > r vz> as 



f(6?r (d+z<) if z e>z< & ^«^< 

L (£<)-(*+*<) if Ze<z< o C«»Ct 



(A. 



This result can also be obtained explicitly by performing 
the integral in Eq. (|A.4bl) . In a similar manner, we find 
that the free energy density J- > can be written in terms 
of the scaling function 

F> =b- {d+z > ) — f > (Rb 1/ ",t < b z <,t > b z >) (A.9) 
V> 

in the regime T > r vz> . Since however £^ <C £,r, m 
this regime, we find that the free energy density of the 
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T ~ \r\ uz < 
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FIG. 4: Generic phase diagrams for systems with multiple 
dynamics. Subfigure (a) and (6) shows the cases of decoupled, 
vz< > 1, and coupled, vz< < 1, multiple dynamic scaling, 
respectively. The regime (II) arises due to the presence of 
multiple thermal scales. 



mode with the larger dynamical exponent always scales 



for T > 



as 



-(d+z>) 



(A.10) 



This can also be derived by explicit integration of the RG 
trajectory (|A.4a|) . 

Thermodynamic observables are defined as derivatives 
of the total free energy density T = T < +F > with respect 
to either the temperature T or the tuning parameter r. 
Here, T and r are the starting values of the RG trajecto- 
ries, which are used in Eqs. (|A.4[) to determine the free 
energy. In the following, we focus on the critical, i.e., the 
most singular contributions to thermodynamics close to 
the quantum critical point. The expression of the flow- 
ing mass r(b) appearing in the RG trajectories for the 
free energy densities can be obtained by integrating the 
corresponding RG equation (|13|) given in the main text. 
For RG stages b 3> &wf ^> 1 with 6wf being the RG 
scale at which the Wilson-Fisher fixed point is reached, 
we obtain that 

r(b) = r (U WF /U f 9 b 1 ' 1 ' (A.ll) 

+ KUwF V< T (b^b^-V'-lfA G(b-b>) , 
\jv — z < \ J 

where Uq = A~ e u /?7< with e = 4 — d — z< being the 
bare reduced interaction and U^j-p its Wilson-Fisher fixed 
point value, and the constant K = A d ~ 2 tt d /(3(2n) d ). A 
derivative with respect to tuning parameter r therefore 
acts on the running r(b) as 



dr(b) 
dr 



= (U WF /U Q ) 4/9 



(A.12) 



and similar terms arise for derivatives with respect to the 
temperature if the RG flow reaches the extended quan- 
tum to classical crossover regime for sufficiently high tem- 
peratures T. 

Each of the two parts, J-> and T < , of the free energy 
density possesses its own crossover between a high and 
low temperature limit. They are determined by the com- 
parison of the correlation length atT = 0, £ r =r _!y , with 
the corresponding thermal length. For J- > , the crossover 



occurs at T ~ r vz> , while it is located at T ~ r VZl - > 
with zr.\ — max{z < ,z e } for J 7 < . These crossovers are 
directly reflected in the behavior of thermodynamic ob- 
servables. We find that the critical part of the specific 
heat, c = —Td\T, is only sensitive to the lower crossover 
line between regime (I) and (II), see Fig. |4j 




if T« r uz > , (I) 

if T » r vz> , (II) + (III) 



(A.13) 



It is always dominated by the mode with the larger dy- 
namical exponent z > because the latter has the larger 
phase space. 

The critical part of the compressibility, n = —d^T, on 
the other hand, is only sensitive to the upper crossover 
line between regime (II) and (III), 



r i/(d+z<)-2 



if T« r uz v , (I) + (II) 



T {d+z < -2/v)/z(,- 1 jf J 1 > r "*(0 , (HI) 



(A.14) 



with Z(.) = max{2<,z e }. The most singular part of the 
compressibility is attributed to the fluctuations with the 
smaller dynamical exponent z < . 

Most interesting is the thermal expansion a = dxd r J-, 
as it is sensitive to both crossovers. We find that the be- 
havior of a at lowest temperatures T < r vz> is dominated 
by the mode with the larger dynamical exponent z > . At 
higher temperatures, the scaling depends on whether de- 
coupled or coupled multiple dynamic scaling is obtained. 
For decoupled scaling present for vz < > 1, the critical 
part reads 



a 



' Tr u{d-i/u- z> ) if T<^r vz >,{\) 
r v{d-i/v) ^ r v Z< r vz> t (n) 

T (d-1M/ Z< jf T^>r uz < , (III). 



(A.15) 

For coupled scaling, z/z< < 1, on the other hand, we 
find that the thermal expansion is sensitive the emergent 
dynamical exponent. We obtain 

' T r u{d-i/v- z> ) jf T<^r vz >,(T) 

jii/(i/z e )-l r i>(d+*<)-2 jf r vze ^> r uz> , (II) 

T (d+z < -z e -i/v)/z e if r»r^,(III). 

(A.16) 

In order to check that the leading behavior indeed 



matches at the lower crossover T 



z> , one needs the 



explicit definitions of v and z e as given in Eqs. (|13j) and 
(fTB)) of the main text. In conclusion, we find that the 
thermal expansion is sensitive to both crossovers in the 
phase diagram of Fig. [4] in agreement with the study of 
Ref . i. 

In the end, we would like to discuss the curvature of 
the upper crossover line between regime (II) and (III). 
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For decoupled and coupled multiple dynamic scaling the 
exponents fulfill the following relations, 



decoupled : v z < > 1 
coupled : v z < < v z e < 1 . 



(A.17) 



As a result, the upper crossover line is concave for cou- 
pled while it convex for decoupled scaling. This serves as 
a convenient necessary condition to identify the type of 
scaling at work, which can be used in experiments. 
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